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Abstract 

Context. White dwarfs are the fossils left by the evolution of low- and intermediate-mass stars, and have very long evolutionary 
timescales. This allows us to use them to explore the properties of old populations, like the Galactic halo. 

Aims. We present a population synthesis study of the luminosity function of halo white dwarfs, aimed at investigating which informa¬ 
tion can be derived from the currently available observed data. 

Methods. We employ an up-to-date population synthesis code based on Monte Carlo techniques, that incorporates the most recent 
and reliable cooling sequences for metal poor progenitors as well as an accurate modeling of the observational biases. 

Results. We find that because the observed sample of halo white dwarfs is restricted to the brightest stars only the hot branch of the 
white dwarf luminosity function can be used for such purposes, and that its shape function is almost insensitive to the most relevant 
inputs, like the adopted cooling sequences, the initial mass function, the density profile of the stellar spheroid, or the adopted fraction 
of unresolved binaries. Moreover, since the cut-off of the observed luminosity has not been yet determined only lower limits to the 
age of the halo population can be placed. 

Conclusions. We conclude that the current observed sample of the halo white dwarf population is still too small to obtain definite 
conclusions about the properties of the stellar halo, and the recently computed white dwarf cooling sequences which incorporate 
residual hydrogen burning should be assessed using metal-poor globular clusters. 
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1. Introduction 


White dwarfs are the evolutionary remnant of stars of interme¬ 
diate and low masses at the zero-age main sequence. The upper 
limit for a main sequence star to evolve to a white dwarf is still 


(iBecker & IbenI 1979 

19801: Mivaii et all 198C 

; Renzini & Voli 

198 It Nomotolll984i 

Garcia-Berro et al.ll 19971 

Poelarends et al. 


200^. Thus, given the shape of the initial mass function it is ex¬ 


pected that the vast majority of the remnants of the evolution of 
single stars will be white dwarfs. Since white dwarfs are numer¬ 
ous, have well studied properties (lAlthaus et al.ll201^ . and have 
long evolutionary timescales, they are the most suitable tool to 
study the properties of old populations, like the Galactic stellar 
spheroid. Moreover, our knowledge of the physics controlling 
the evolution of white dwarfs relies on solid grounds, since the 
basic principle of their evolution is a well understood and rela¬ 
tively simple cooling process. Although this basic principle of 
the theory of white dwarf cooling has remained unaltered dur¬ 
ing the last decades, we now have very sophisticated and accu¬ 
rate stellar evolutionary models that allow us to perform precise 
cosmochronology, and to characterize the ensemble properties 
of several white dwarf populatio ns, like those of the Galactic 
disk - see ICoiocaru et alJ (l2014l) . and references therein, for a 
recent work on this subject - and of the system of Galactic open 
(iGarefa-Berro et ani201(ll: Bellini_et^ 201(1 Bedin_et^ 201T) 


and globular clusters (iHansen et al.F 20021 : iHansen et all 2013 

iGarefa-Berro et alJl2014ll. 


The population of white dwarfs in the Galactic stellar halo 
has been the subject of an increase d interest since the first 
observationa l and t heoretical studies (iMochkovitch et aklllQQOl : 
iLiebert et al.l 1198^ . Perhaps, one of the most important rea¬ 
sons of this interest in halo white dwarfs is their possible con¬ 
tribution to the dark matter conten t of our Galaxy - see, for 
inst ance, Qppenheimer et all (1200 ll) for an observational work, 
and iTorres et al.l ( 20021) for a theoretical study. However, due to 


the very low space densities and the intrinsic faintness of the 
population of white dwarfs of the Galactic spheroid, their detec¬ 
tion has proven to be a difficult endeavour. Moreover, opposite 
to what occurs with main sequence stars which can be classified 
according to their metallicity, the atmospheres of white dwarfs 
are devoided of metals. This is due to their high surface gravi¬ 
ties and long evolutionary timescales, which allow gravitational 
diffusion to be very efficient at settling the metals resulting from 
the previous evolutionary history at the base of the partially de¬ 
generate envelope. All these physical processes make halo white 
dwarfs indistinguishable from disk ones. Hence, the only obser¬ 
vational method to detect white dwarfs belonging to the Galactic 
spheroid not hampered by relevant technical difficulties relies 
on identifying them on the basis of large proper motions, as ra¬ 
dial velocities cannot be determined accurately due, again, to 
the large surface gravities, which translates into a sizable gravi¬ 
tational redshift of the spectral features that cannot be neglected 
and is difficult to measure. Additionally, the absence of spec¬ 
tral lines at the very low luminosities of the coldest halo white 
dwarfs also prevents an accurate characterization of the faintest 
population of halo white dwarfs. All this, in turn, reduces consid- 
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erably the size of the observational sample, since at present large 
volumes cannot be probed, and we are limited to study nearby 
halo white dwarfs. 

Nevertheless, recent observational attempts to empirically 
determine the luminosity function of halo white dwarfs have 
come to a succes s, and we now have a reliable sample o f 
halo white dwarfs (iHarris et al.ll2006HRowell & Hamblvll201 ih . 
to which the theoretical works can be compared. Comparing 
the results of the theoretical models with the available ob¬ 
served sample of halo white dwarfs is an important task, as 
we now also have accurate white dwarf cooling tracks for 
whit e dwarfs descending from very low-metallicity progeni¬ 
tors (iMiller Bertolami et al.ll2013t lAlthaus et al.ll2015l) that im- 
pr ove upon those us ed in the early and pion e ering calculations 
of llsern et al.l (1199^ and Garcfa-Berro et al. (l2004l) . and in the 
more recent one of Ivan Oirschot et alJ (12014 ). These evolution¬ 
ary sequences have been evolved self-consistently from the zero 
age main sequence, through the red giant and thermally pulsing 
AGB phases to the white dwarf regime, and have revealed the 
important role of residual hydrogen burning in the atmospheres 
of low-mass white dwarfs, a physical process that needs verifica¬ 
tion. Finally, these kind of works are also of crucial importance 
to pave the road to future studies of the large population of halo 
white dwarfs that the European astrometr ic mission Gaia will 
unveil in the next years (iTorres et al.ll2005l) . 

Our paper aims at producing synthetic samples of the pop¬ 
ulation of halo white dwarfs using the most up-to-date physical 
inputs and prescriptions for the Galactic spheroid and compare 
them with the current observational data. It is organized as fol¬ 
lows. In Sect. 12.11 we briefly describe the numerical tools em¬ 
ployed in this work. It is followed by Sect. [2 where we first dis¬ 
cuss the effects of residual hydrogen burning, of the adopted ini¬ 
tial mass function, of the assumed density profile for the Galactic 
halo, of a population of unresolved binary white dwarfs, and 
thosee of the star formation history. Finally, in Sect. |4] we sum¬ 
marize our calculations and we draw our conclusions. 


2. The population synthesis code 

2.1. A brief description of the numericai set-up 


As in ou r prev i ous works (Garcfa-Berro et al.] 199^ 


Torres et ^ 2002: Garcfa-Berro et al.r ~ 2004 : [Torres et dl 


2005t ICoiocaru et Sll2014ll . we use a Monte Carlo population 


synthesis code, in this case adapted to model the halo popula¬ 
tion. In the following we describe the most important inputs of 
our standard model. 

We initially produced a large number of synthetic main se¬ 
quence stars, located according to a isothermal sphere density 
model: 


p(r) oc 


H- Rq 
H- 


( 1 ) 


where a ^ 5 kpc is the core radius, and Rq = 8.5 kpc is the 
galactocentric distance of the Sun. We assigned to each syn¬ 
thetic star a value for the mass at the zero-age main s equence 
random ly generated from the initial mass function of iKroup^ 
(1200 ill , and a time of birth, randomly assigned within a burst of 
constant star formation lasting for 1 Gyr. The velocities of simu- 
lated halo stars were rando mly drawn from normal distributions 
dBinnev & Tremai'n9ll987l) : 
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where cTr and cr, - the radial and the tangential velocity disper¬ 
sions, respectively - are related by the following expression: 
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which, to a first approximation, leads to cry - o-f - Vd V2. 
The velocity dispersions cTr and cr , are those determined by 
iMarkovic & Sommer-Lars^ (Il997l) . For the calculations re¬ 
ported here we adopted a circular velocity 14 = 220 km/s. From 
these velocities we obtained the heliocentric velocities of each 
simulated star by adding the velocity of the local standard of rest 
(FSR) vlsr = -220 km/s, and the peculiar velocity of the Sun. 

Next we computed the main sequence lifetime for each pro¬ 
genitor star, adopting a set of evolutionary sequences with metal- 
licity Z = 0.0001, which together with the age of the population 
(for which in our reference model we adopted 14 Gyr), and the 
progenitor mass, allowed us to determine which stars have had 
time to become white dwarfs at present time. We then obtained 
the corresponding masses and cooling ages for each simulated 
white dwarf. It is worth mentioning here that all the evolution¬ 
ary sequences - that is, the progenitor evolutionary sequences 
and the white dwarf cooling ones - employed in our work are 
those of lAlthaus et al.l (12015h . which were obtained from fully 
ev olutionary calculations , and e xpand the previous calculations 
of iMiller Bertolami et akl (1201 3h . Hence, the main sequence life¬ 
times, the relationship linking the progenitor and the white dwarf 
masses, and the cooling ages are all self-consistently computed 
using an homogenous evolutionary framework. This represents a 
clear improvement over the most recent calculations of this kind, 
as we employed self-consistent evolutionary models of the right 
metallicity that incorporate state-of-the-art prescriptions for all 
the relevant physical processes. Our calculations incorporate a 
fraction of 20% of non-DA white dwarfs, for which we employ 
theoretical cooling sequences for white dwarfs with pure helium 
atmospheres. We elaborate on the cooling tracks employed here 
in Sect. 12.21 Using these values we derived the stellar parame¬ 
ters of each white dwarf in the synthetic sample. Namely, we 
computed its luminosity, effective temperature, surface gravity 
and magnitudes in the different pa ssbands. A standard model of 
Galactic absorption was also used (iHakkila et al.ll 19971) to obtain 
reliable apparent magnitudes. 

Our synthetic white dwarf sample is then passed through a 
series of Alters which mimic the selection criteria employed to 
observationally select halo white dwarfs in a real sample. These 
Alters are described in detail in Sect. 12.31 After this procedure 
is followed the white dwarf luminosity function can be com¬ 
puted, except for a normalization factor. We chose to normal¬ 
ize the theoretical results to the density of white dwarfs in the 
highest density bin with finite error bars of the observational lu¬ 
minosity function, Mboi = 15.75. This is, in fact, equivalent to 
normalize the luminosity function to the total population den¬ 
sity, given that this bin practically dominates the stellar counts. 
We remark at this point that in our fiducial model only single 
white dwarfs were considered, however we also explored mod¬ 
els with a fraction of unresolved binaries in our calculations. It is 
also worth mentioning here that our simulations include as well 
a careful exploration of the effects of other inputs, which will be 
further explained in Sect. [2 


2.2. Cooiing tracks 

White dwarf progenitors in the Galactic halo are characterized 
by a significantly low metallicity. In the Solar vicinity, the halo 
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Table 1. Number of synthetic white dwarfs that survive the dif¬ 
ferent observational cuts for a typical Monte Carlo realization of 
our standard model. 


Filter 

Nwd 

% 

Initial sample 

592 199 

100 

/^min cut 

8 952 

1.5 

12 < r59p < 19.75 

111 

0.02 

RPMD cut 

111 

0.02 

Uan > 200 km 

77 

0.01 


metallici ty distribution function pe aks at rFe/H1~ -1. 5 dex. 
Actually. iFrebel & Norrisl(l2013h and ICarollo et all (1201 Oh found 
that the Galactic halo has a dual population. The first of these 

halo populations peaks at [Fe/H]-1.6 dex, whereas the second 

one peaks at [Fe/H]-2.2 dex. All in all, it is clear that to ade¬ 

quately capture the essential properties of this metal-poor popu¬ 
lation, a set of cooling sequences of white dwarfs with hydrogen- 
rich atmospheres descending from low-metallicity progenitors is 
needed. 

We interpolate th e cooling tim e s usin g the set of evolu¬ 
tionary sequences of lAlthaus et al.l (l2015h . These cooling se¬ 
quences have been computed considering stable, residual hydro¬ 
gen shell burning in white dwarf atmospheres during the white 
dwarf stage, although they also provide a set of cooling tracks 
in which this physic al mechanism is disreg a rded. T his is an im¬ 
portant issue, since [Miller Bertolami et al.l (1201 3h showed that 
although in most cases residual hydrogen burning is not a sig¬ 
nificant source of energy, for white dwarfs with hydrogen at¬ 
mospheres descending from progenitors with very low metal¬ 
licity it can become a dominant source of energy, and can de¬ 
lay significantly white dwarf cooling. This effect is more notice¬ 
able for low-mass white dwarfs with luminosities ranging from 
log(L/LQ) = -2 to -4. As mentioned, we consider that the adop¬ 
tion of this set of sequences represents a clear improvement with 
respect to the most recent calculation of the luminosity of halo 
white dwarfs (Ivan Oirschot et al1l2014ll . which employed evolu¬ 
tionary sequences for progenitors of Solar metallicity. 

Although the evolutionary sequences for white dwarfs with 
hydrogen-rich atmospheres adopted in this work are a clear im¬ 
provement over previous attempts to model the population of 
single white dwarfs in the Galactic halo, a cautionary remark 
is in order here. There is solid evidence t hat old stellar sys¬ 
tems exhibit an enhance ment of a elements (lAller & Greensteinl 
ll96Ct IWallersteinlll962h . While such an enhancement has yrr- 
tually no effects on the eyolutionary timescales of initially low 
mass stars, they can play a role in the eyolution of intermediate 
mass stars. In particular, the resulting total metallicity is larger 
than the one obtained by assuming a Solar-scaled composition 
and, due to the increase of the oxygen abundance, the global 
abundance of carbon, nitrogen and oxygen is larger than the cor¬ 
responding Solar-scaled abundance. This, in turn, has an effect 
on the eyolutionary timescales of the progenitor stars of typi¬ 
cal white dwarfs. A more rigorous treatment of low-metallicity 
stars should require the inclusion of a-enhanced initial chemi¬ 
cal compositions to compute the stellar sequences. Our eyolu¬ 
tionary sequences do not take into account this enhancement, 
but we estimate that the effects of including it in the calcula¬ 
tion of the white dwarf luminosity function is limited. In partic¬ 
ular, for the metallicities and progenitor masses releyant for our 
study, we haye checked that the differences of the progenitor life¬ 
times and of the resulting white dwarf masses between the yal- 


ues obtained with solar-scaled metallicities and those obtained 
with an a-enhanced che mical composition are smaller than 1 % 
(iPietrinferni et al.ll2006h . Hence, our results are almost insensi- 
tiye to the adopted metal ratios. 

Finally, for more massiye oxyg en-neon white dwarf s we 
employed the coolin g sequences of lAlthaus et al.l (l2005l) and 
lAlthaus et alJ (l2007h . whereas for white dwarfs with pure he¬ 
lium a tmospheres we used the cooling tracks of [Bergeron et al.l 
(1201 Ih . In both cases the white dwarf eyolutionary sequences 
correspond to progenitors of Solar metallicity. This, of course, is 
not a self-consitent treatment, but neyertheless we judge that the 
effects on the computed white dwarf luminosity functions should 
be modest, see below for a detailed discussion. 

2.3. The observational sample and its observational cuts 

We compare our results with the most recent and statistically 
releyant obseryational h alo white dwarf luminosity function 
(iRowell & Hambl^l201 Ih . This obseryational luminosity func¬ 
tion was deriyed from a sample of 93 halo white dwarfs de¬ 
tected in the SuperCosmos Sky Suryey (SSS). The SSS is an 
adyanced photographic plate-digitizing machine, using plates 
taken with the UK Schmidt telescope (UKST), the ESO Schmidt 
telescope, and the Palomar Schmidt telescope. The resulting cat¬ 
alogs were compiled by digitizing seyeral generations of pho¬ 
tographic Schmidt plates. The suryey uses a ph otometric sys¬ 
tem t hat has three passbands: bj, and In (iHambly et al.l 
20011). Employing data from seyeral generations of plates, 
Rowell & Haniblyl (1201 ll) constructed a catalog of ~ 10,000 
white dwarfs with magnitudes down to rsgp ~ 19.75, and proper 
motions as low as p ~ 0.05 yr ', coyering nearly three quar¬ 
ters of the sky. Using strict yelocity cuts, the authors isolated 
subsamples of white dwarfs belonging to the thin disk, the thick 
disk and the halo populations, and presented obseryational white 
dwarf luminosity functions for each one of these populations. 

In our study we distinguish between the complete sample 
of synthetic halo white dwarfs and a restricted sample. The lat¬ 
ter is obtained by replicating the obseryational selection crite¬ 
ria adopted to deriye the obseryed halo white dwarf luminosity 
function of the SSS. Eirst, a proper motion cut depending on the 
bj magnitude is applied. This proper motion cut is giyen by the 
following expression; p > 5((T™’‘(hj) -i- 0.002), where cr^ is the 
standard deyiation in the proper motion measurements. Also a 
magnitude cut is imposed, 12 < r^gp < 19.75. Next, a cut in 
the reduced proper motion diagram is performed, selecting only 
objects below and blueward of a reduced proper motion corre¬ 
sponding to Vtan = 30 km s“'. Lastly, in order to separate the 
halo population, a tangential yelocity cut is used. Specifically, 
we only select stars with tangential yelocities Vt^n > 200 km s“*. 
Einally, we also impose an upper limit on the tangential yeloc¬ 
ity of 400 km s“' to preyent selecting stars with yelocities larger 
than the escape yelocity of the Galaxy. 


3. Results 

In this section we compare t he results of our simulatio ns to the 
halo luminosity function of iRowell & Haml^ (1201 ll) . and we 
study the sensitiyity of the theoretical white dwarf luminosity 
function to different model inputs. 

To start with, we discuss how the obseryational selection cri¬ 
teria affect the size of the synthetic samples. This is done with 
the help of Table [1] In this table we list for our reference model 
the number of white dwarfs of the original synthetic sample (first 
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Figure 1. Halo white dwarf luminosity function for our fidu¬ 
cial Galactic model. The top panel shows the theoretical white 
dwarf luminosity function obtained when the cooling sequences 
incorporating residual hydrogen burning are employed (open 
ci rcles). We also show the observed halo luminosity function 
of iRowell & Hambhl (1201 ih - solid squares. The bottom panel 
shows the residuals between the observed luminosity function, 
and the theoretical calculations, A = 2(A^obs-A^sim)/(A^obs+A^sim)- 
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Figure 2. White dwarf luminosity functions when only 
hydrogen-rich synthetic white dwarfs are generated in the Monte 
Carlo simulation. The bottom panel shows the residuals between 
our standard model, in which a fraction of 20% of non-DA white 
dwarfs was employed, and that in which this percentage is zero, 
A = 2(Astd-A^no-DA)/(A^std+A^no-DA) — solid squares and hollow 
circles, respectively. 


row) and in subsequent rows we list the number of white dwarfs 
that survive the different cuts. As can be seen only 1.5% of the 
synthetic stars survive the proper motion cut. After applying the 
magnitude cut we are left with 111 synthetic stars, representing 
about 0.02% of the original sample. For this particular realiza¬ 
tion the reduced proper motion cut does not decrease further the 
number of simulated white dwarfs, whereas the filter in tangen¬ 
tial velocities reduces even more the number of simulated stars 
to about 0.01% of the original sample, to 77 white dwarfs, a 
number comparable with that observationall y found. _ 

To compare our results with those of Ivan Oirschot et al] 
(l2014l) we only culled white dwarfs using the tangential veloc¬ 
ity cut, as they did. Using only this selection criterion the size 
of the restricted sample is much larger. In particular, when this 
procedure is adopted it results in a restricted sample which is 
63% of the initial one. Obviously, the advantage of such a large 
synthetic sample is that it is comparable to the complete one, 
producing a smooth luminosity function that faithfully preserves 
the intricacies of the adopted model. However, the main draw¬ 
back of adopting just this selection criterion is that the resulting 
sample is ultimately not directly comparable to the observational 
one, as is ours. 

The top panel of Fig. [T] shows the white dwarf luminosity 
function of our reference mo del - open circles - and the ob¬ 
served luminosity function of iRowell & Hambhl (1201 ll) - solid 


squares - while in the bottom panel the corresponding residu¬ 
als are displayed. As can be seen, the agreement between the 
theoretical results and the observed data is very good. Our fidu¬ 
cial model reproduces not only the observed slope of the white 
dwarf luminosity function, but also accounts for the scarcity of 
halo white dwarfs at very low luminosities (Mboi > 17). This 
indicates that our Monte Ca rlo code correctly reproduc es the se¬ 
lection criteria employed bv IRowell & Hamblvl (1201 lit . 

In a second step, we have checked the sensitivity of our syn¬ 
thetic white dwarf luminosity function to our choice of cool¬ 
ing sequences for massive, ONe white dwarfs, and non-DA 
white dwarfs. We recall that for these stars we employed a set 
of cooling sequences of Solar metallicity. Specifically, we as¬ 
sessed which is the final number of these white dwarfs in the re¬ 
stricted sample - that is, once the observational selection criteria 
are taken into account - and we found that in a typical Monte 
Carlo realization only one of these white dwarfs, at most, sur¬ 
vives the successive selection cuts. The most stringent observa¬ 
tional cut is the magnitude cut, rsgp ~ 19.75. In most Monte 
Carlo realizations none of these white dwarfs survives this cut. 
Additionally, we mention that even if this cut is not employed, 
the proper motion cut eliminates from the final sample almost 
99.5% of ONe white dwarfs. Thus, there are very few ONe 
white dwarfs in the final sample. The reason for this behav¬ 
ior is twofold. First, these white dwarfs are very scarce, since 
their formation is strongly inhibited by the shape of the initial 
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Figure 3. Same as Fig.[T]for the case in which residual hydrogen 
burning is not considered. The bottom panel shows the residuals 
between the luminosity function computed using our standard 
cooling sequences and that obtained when residual hydrogen 
burning is artifically ignored, A = 2(A^std-A^no-H)/(A^std+A^no^H)- 


Figure 4. White dwarf luminosity functions when different ini¬ 
tial mass functions are considered for the stellar spheroid. The 
bottom panel shows the resi duals between o ur standard model 
and those obtained wh en the Salpetei] ( 1955h and the top-heavy 
initial mass function of ISuda et al. ( 2013h are used, A = 2(A^std - 
NiMF)/{Ns[d + Nimf) — void circles and triangles, respectively. 


mass function. Thus, not surpringsily, they contribute little to 
the white dwarf luminosity function. The second reason is that 
being these white dwarfs made of oxygen and neon, their heat 
capacity is sm aller than that of a carbon -oxygen white dwarf of 
the same mass (iGarcia-Berro et af .Il997h . and consequently cool 
faster. Accordingly, these white dwarfs contribute essentially to 
the faintest bins of the luminosity function, which is not probed 
by observations, because it is excluded by the magnitude cut. 
In summary, we conclude that the influence of adopting a set of 
cooling sequences of Solar metallicity for ONe white dwarfs is 
negligible. 

To assess the influence of adopting a set of cooling sequences 
of Solar metallicity for non-DA white dwarfs we ran an ad¬ 
ditional simulation in which the percentage of non-DA white 
dwarfs was set to zero, and consequently all the synthetic stars 
had hydrogen-rich atmospheres. We then computed the residu¬ 
als between the resulting white dwarf luminosity function and 
that obtained with our reference model, for which this percent¬ 
age is 20%. The results are shown in Fig. |2] As can be seen 
the differences are small, although not negligible. As a matter 
of fact, the space density of hot white dwarfs is smaller in the 
case in which only synthetic DA white dwarfs are generated. 
However, this is a consequence of the normalization procedure. 
Hydrogen-defficient white dwarfs have cooling sequences that 
resemble those of a black body, whereas the atmospheres of DA 
white dwarfs are more transparent. Consequently, at low tem¬ 
peratures non-DA white dwarfs cool faster than DA ones, and 
the percentage of non-DA white dwarfs increases for decreasing 


luminosities, so these white dwarfs accumulate at luminosities 
close to that of the peak of the theoretical luminosity function 
and even smaller. However, the number counts of white dwarfs 
in the luminosity bins close to the peak of the luminosity func¬ 
tion dominate the total number counts of white dwarfs in the 
synthetic sample. Thus, since the total number of white dwarfs 
in any Monte Carlo realization must be kept constant and, more¬ 
over, must be close to the observed value, the hot branch of the 
luminosity function is depleted in the case in which non-DA 
white dwarfs are not generated. Nevertheless, we emphasize that 
because white dwarf cooling sequences of low metallicity for 
non-DA white dwarfs are not available, it is clear that this proce¬ 
dure largely overestimates the impact of adopting a set of cool¬ 
ing sequences of unappropriate metallicity. Thus, we conclude 
that the possible effect of adopting a set of cooling sequences of 
Solar metallicity for non-DA white dwarfs is limited. 

Next, we assess the sensitivity of these results to the most 
relevant inputs of our model. In particular, we first discuss if 
the adopted cooling tracks for carbon-oxygen white dwarfs with 
hydrogen-rich atmospheres could change this picture. In a sec¬ 
ond step, we study whether a different choice of the adopted ini¬ 
tial mass function could affect our results. Later, we evaluate if 
a different halo model could have a noticeable influence in our 
calculations. Finally, we also study if different percentages of 
unresolved binaries vary the shape of the white dwarf luminos¬ 
ity function. We end our assessment comparing our theoretical 
resuls for different ages of the stellar halo. 
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Figure 5. White dwarf luminosity functions for different density 
profiles of the stellar halo. The bottom panel shows the residuals 
between our standard mo del and those obtained w hen the density 
profile of a triaxial halo (lHelmill2004t) and that of iNavarro et alJ 
(1 19961) are employed, A = 2(A^std-A^p)/(A^std+A^p) — void circles 
and triangles, respectively. 
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Figure 6. Same as Fig.|5]but for different fractions of unresolved 
binaries. The bottom panel shows the residuals between our stan¬ 
dard model and those obtained when the adopted fractions of un¬ 
resolved binaries are 20% and 40%, A = 2(A^std - A^bin)/(A^std + 
A^bin) — void circles and triangles, respectively. 


3.1. Hydrogen burning 

It has been shown (iMiller Bertolami et alj201^ that residual hy¬ 
drogen burning can have a significant impact on the cooling pro¬ 
cess of white dwarfs with progenitors of very low metallicity, the 
effect being more noticeable for low-mass white dwarfs (those 
with masses between 0.5 and 0.6 Mq). Since low-mass white 
dwarfs contribute to all the luminosity bins of the hot branch 
of the luminosity function, and since the shape of the luminos¬ 
ity function is directly related to the cooling rate, it is natural 
to ask ourselves whether a different choice of cooling sequences 
could affect the slope at moderately high luminosities. We check 
this using the two d ifferent sets of cooling tracks described in 
lAlthaus et al.l (l2015h . The first of these sets is the one used in our 
reference model, and considers residual nuclear burning, while 
the second one does not take into account nuclear reactions (as 
it occurs for white dwarf progenitors with Z > 0.001). 

In Fig. |3] we present the resulting white dwarf luminosity 
function for the halo population when the cooling sequences 
in which residual hydrogen burning is artificially ignored are 
adopted. This luminosity function should be compared with that 
shown in Fig. [T] As can be seen, the only difference between 
both sets of theoretical calculations is that for the case in which 
the cooling sequences incorporating residual hydrogen burning 
are employed there is a small plateau between Mboi = 12 and 
14, which is absent in the case in which no residual hydrogen 
burning is considered. This plateau reflects the slow-down of 


cooling due to the release of energy of residual hydrogen burn¬ 
ing. However the differences between both calculations are mi¬ 
nor, and the currently available observational luminosity func¬ 
tion (which is derived using only ~ 100 white dwarfs) does not 
allow to draw definite conclusions about the real existence of 
residual nuclear burning. 


3.2. Initial mass function 


As mentioned, we also test the influence that the adopted ini¬ 
tial mass function may have in our results. Since the formation 
timescale of the stellar halo is short, it is straightforward to show 
that when a burst of negligible duration is adopted the luminosity 
function is given by 


N(L) oc 


dn 

dM\,o\ 


dn dM dM 

-oc ^{M) - 

dMdMboi dM^x 


(4) 


In this expression n stands for the space density, and O for the 
initial mass function (IMF). Thus, it is clear that the adopted 
initial mass function should influence the shape of the luminosity 
function. 

In order to test the influence of the IMF on the luminos¬ 
ity function we employ three parametrizations. The first one is 
that used in our f iducial model, namely the so-called “universal” 
mass function of iKroun^ (1200 ll) . For the mass ranges relevant to 
our study this IMF is totally equivalent to a two-branch power 
law with exponent -a, with a = 1.3 for 0.08 < M/Mq < 0.5 
and a - 2.3 for M/Mq > 0.5. We also compute theoretical 
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Figure?. Same as Fig. |5] but for different durations of the ini¬ 
tial burst of star formation. The bottom panel shows the resid¬ 
uals between our standard model and those obtained when the 
adopted durations of the initial burst of star formation are 0.5 
and 2.0 Gyr, A = 2(A^std - A^A/)/(A^std + N^t) — void circles and 
triangles, respectively. 


white dwarf lum inosity functions adopting the classical IMF of 
ISalneteii (Il955h . which is a power law with index a - 2.35. 
Finally we also adopt a top-heavy IMF; 


(1)(M) = — exp 


-log(M/;u) \ 
2cr2 I 


(5) 


In this expr ession u - I OMiT) and cr — 0.44. This IMF was in¬ 
troduced bv ISuda et al.l d2013h . and is dominated by high mass 
stars. It has been found that this IMF better reproduces the 
characteristics of metal-poor populations, namely those with 
[Fe/H]< -2. 

The corresponding luminosity functions for these IMFs are 
shown in the top panel of Fig. [H and their respective residuals 
with respect to our fiducial model are shown in the bottom panel 
of this figure. As can be seen, there are no notice able diff e rences 
between t he calculations i n which the IMF of iKroupal (1200 Ih 
and that of lSalneteil (Il955h are employed. The reason for this is 
that in the relevant luminosity range the slope of both IMFs is 
very similar. We note, however, that when the top-heavy IMF of 
ISuda et ahl (1201 3h is used the luminosity function presents a drop 
in the space density at large luminosities. This deficit of bright 
white dwarfs is quite apparent, but it is marginally consistent 
with the observed data. 


3.3. Density profiles 

Another possible concern would be the adopted density profile 
for the stellar halo. As explained in Sect. 12.11 in our reference 
model we adopted the density profile of the classical isother¬ 
mal sphere, but there are other density profiles that are worth 
studying. Accordingly here we study how this choice affects our 
results. To do this we first adopted a triaxial oblat e halo mode l, 
which is based in a logarithmic dark halo potential (lHelmil2004l) . 


V^^vo^ln(R^+z^/q^ + d^) 

being the corresponding density distribution: 




' vg ^ jlq^ + \)cfi +R^ + {2- q-^)z^ 
AnGq^ ^ {cP- +R'^ + z^q^^P 


( 6 ) 


(7) 


In this expression we have adopted d - 12 kpc and vq = 
131.5 km/s, which gives a circular velocity of the Sun of 
229 km/s, and an oblatness parameter q = 0-8. T he third and last 
profile we used is that of ((Navarro et al.lll99^ . a widely used 
one; 



with rj = 18 kpc. 

As Fig. |5] reveals, the differences between the luminosity 
functions computed using these three different density profiles 
for the stellar halo are totally negligible. This is because the sam¬ 
ple of halo white dwarfs of iRowell & Hambl^ (1201 Ih is local, 
whereas the differences between the three model profiles should 
be prominent at large distances. 

3.4. Unresolved binaries 

One of the potential problems when calculating the observed lu¬ 
minosity function for single stars are unresolved binary white 
dwarfs, since they compute as single stars, and hence can result 
in a different shape of the luminosity function. T his has been 
proven to be the case in some Galactic clusters (iBedin et al.l 
120081 : iGarcfa-Berro et all l2010l) . It is therefore interesting to 
check the effect that a certain fraction of unresolved binaries 
can have on the theoretical luminosity function. In order to test 
this we compute a new set of simulations, based on our fidu¬ 
cial model, and increasing the fraction of unresolved binaries. 
We remind that in our reference model no unresolved binaries 
were considered. As for the distribution of secondary masses we 
adopted a model in which the masses of both components are 
uncorrelated. 

Fig. 12 shows the result of this numerical experiment when 
the fractions of unresolved binaries are, respectively, 20% and 
40%. As can be observed in this figure, increasing the fraction 
of unresolved binaries considered in the sample does not result in 
any noticeable change, but in a slight reduction of the number of 
white dwarfs populating the brightest luminosity bins. The rea¬ 
son for this can easily be explained. Since low-luminosity white 
dwarfs have longer evolutionary timescales the low-luminosity 
bins have also large space densities. Consequently, unresolved 
binaries also concentrate in the luminosity bins with the largest 
densities, and thus the bright luminosity bins are less populated. 
Since we normalize our theoretical luminosity function to the 
observed luminosity bin at Mboi = 15.75, the result is that the 
bright branch of the theoretical luminosity function is depleted. 
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Figures. Same as Fig.|5]but for four various merger episodes, of two strengths at two different times. The bottom panels show the 
residuals between our standard model and those obtained when the impact of a merger episode is analyzed, A = 2(A^std-A^mer)/(A^std+ 
A^tner)- See text for details. 


Nevertheless, the differences are minor even when an unrealis¬ 
tic percentage of 40% of the objects in the synthetic sample are 
unresolved binaries. 


3.5. The star formation history 


Another point of concern is the adopted star formation history. 
This also may have potential effects on the morphology of the 
hot branch of the halo white dwarf luminosity function. To start 
with we discuss the effects of the duration of the initial burst of 
star formation. This is done with the help of Fig. |7] where we 
show the theoretical white dwarf luminosity functions for two 
burst of durations 0.5 and 2.0 Gyr, and compare them with our 
reference model, for which we recall we employed a burst of 
duration 1.0 Gyr. This figure clearly shows that, except for the 
smaller space densities at moderately high luminosities, the dif¬ 
ferences between these two luminosity functions and our fiducial 
one are marginal. Consequently, current observations do not al¬ 
low to discern between different durations of the initial burst of 
star formation. 


Furthermore, a consensus about the origin of the stel¬ 
lar spheroid has not been reached yet. The two main com- 
petin g scenarios - mo nolithic collpase of the protogalactic 


-1 

gas (lEggen et al.l I196 2[). or fo rmation through several merger 
episodes (ISearle & Zinnll97^ - need still to be confronted with 
observations. Hence, it is natural to wonder if these two scenar¬ 
ios leave observable imprints in the shape of the hot branch of the 
white dwarf luminosity function of single halo white dwarfs. To 
this end we conducted an additional set of simulations in which, 
in addition to the initial burst of star formation, we modeled the 


luminosity function when a second burst of star formation oc¬ 
curring some time ago is adopted. Specifically, we ran four ad¬ 
ditional simulations in which a secondary burst of star formation 
occurs at times 4 and 8 Gyr respectively, varying the strength of 
this secondary burst. The metallicities of the secondary bursts of 
star formation were the same adopted for the initial one. This 
choice minimizes the effects of such merger episodes, but the 
effects of the different metallicity of the secondary bursts are ex¬ 
pected to be minor. Specifically, the secondary burst was given 
amplitudes 20% and 40% of the initial one. In all cases, the du¬ 
rations of all the bursts (that is, both the initial and the secondary 
ones) were kept fixed and equal to 0.1 Gyr, while we recall that in 
the standard model a duration of 1 Gyr was adopted. The results 
of this numerical experiment are displayed in Fig. [8] In the left 
panels of this figure we show the results when a secondary burst 
with an amplitude 20% of the initial one is adopted, whereas in 
the right panels the results when the amplitude of the secondary 
burst is increased to 40% of the primary one are displayed. As 
can be seen, the differences are again very small. Thus, unfortu¬ 
nately, the current observational database of halo white dwarfs 
does not allow to distinguish the two aforementioned formation 
scenarios of the stellar halo. 


3.6. Age of the popuiation 

Finally, we ran a set of simulations in which we varied the age 
of the halo population, from 11 to 13 Gyr, and we compared the 
results of these calculations with that obtained in our reference 
model, for which we remind we adopted an age of 13.7 Gyr. We 
show the results of these calculations in Fig.|2 As expected, the 
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Figure 9. Same as Fig.|5]but for different ages of the halo popu¬ 
lation. The bottom panel shows the residuals between our stan¬ 
dard model and the luminosity functions obtained when the age 
of the stellar halo is varied, A = 2(A^std - A^r)/(A^std + Nt)- 


bright branch of the white dwarf luminosity function does not 
depend appreciably on the adopted age of the stellar spheroid. 
Moreover, since the observed luminosity function does not show 
a cut-off the age of the halo population cannot be yet computed 
using the termination of the cooling sequence of halo white 
dwarfs. This is a consequence of the cuts used to select the ob¬ 
served sample, and specifically is caused by the cut in bolomet- 
ric magnitude. The only quantitative assessment about the age 
of the halo that can be made with the available observed data 
is to place a lower limit. This can be done in a simple way, by 
imposing that the dimmest populated luminosity bin of the the¬ 
oretical white dwarf luminosity function is that observationally 
found, at Mboi = 17.25. Using this procedure we find that, al¬ 
though it is not possible to fit the halo age, a lower limit for its 
age of 12.5 Gyr can be safely established. 


4. Conclusions 

In this paper we have revisited the luminosity function of halo 
white dwarfs, in the light of the recently computed white dwarf 
cooling seque nces for low-metallicity progenitors. These cool- 
ing sequences (iMiller Bertolami et alj201^lAlthaus et al.l2015l) 
have been derived evolving self-consistently their progenitors 
from the zero age main sequence, through the red giant and 
thermally pulsing AGB phases to the white dwarf regime, and 
have unveiled the role of residual hydrogen burning in the atmo¬ 
spheres of low-mass white dwarfs. In this sense, it is important to 
realize that these evolutionary calculations superseed those used 
in the early and pioneering calculations of the halo white dwarf 


lumino sity function of llsern et al.l (Il998h and Garcia-Berro^^l.] 
(l2004ll . and in the recent work of van Oirschot et al.l ( 20141) . 
Moreover, in pursuing this endeavour we have employed a 
state-of-the-art numerical code, which incorporates not only the 
most recent advances that allow an accurate description of the 
Galactic halo, but also a detailed implementation of the obser¬ 
vational biases and restrictions, an issue that most theoretical 
calculations do not take into account, thus impeding a sought 
comparison with the observed sample. This is an important is¬ 
sue, as the observed sample of white dwarfs belonging to the 
Galactic spheroid suffers from small statistics. Moreover, given 
that the density of halo white dwars is low and that this popula¬ 
tion is old, hence intrinsically faint, the detection of halo white 
dwarfs is hampered by observational difficulties. Consequently, 
the selection biases are important, and we are restricted to com¬ 
pare the theoretical results of our results with an observational 
sample plagued with uncertainties. 


Since residual hydrogen burning occurs at moderately low 
luminosities, say from log(L/LQ) = -2 to -4, the halo lu¬ 
minosity function could eventually offer an unique possibility 
to test the reliability of these recent cooling sequences. This 
could have important consequences for our understanding of 
how white dwarfs are formed and how their progenitor stars 
evolve in low-metallicity environments, and more specifically 
it could shed light on the occurrence of the third dredge-up 
for metallicities 10“^. We have found that, unfortunately, the 
scarcity of halo white dwarfs at the luminosities at which resid¬ 
ual hydrogen burning occurs prevents us from making a mean¬ 
ingful comparison between the sequences which incorporate this 
physical ingredient and those which do not. Thus, this effort will 
have to wait until we have larger and more reliable samples. 
Alternatively, this could be done using the white dwarf luminos¬ 
ity functions of Galactic globular clusters, of which NGC 6397 
is, perhaps, the leading example. 


Additionally, we have investigated whether or not the ob¬ 
served luminosity function of single white dwarfs can be even¬ 
tually used to learn more about the stellar population of the 
Galactic halo. In particular, we have studied if it can be used to 
constrain the initial mass function of this population, its star for¬ 
mation history and age, to probe different halo density profiles, 
or possibly to discern the fraction of unresolved binaries that 
may contaminate observations. Unfortunately, our calculations 
show that the hot branch of the luminosity function is almost in¬ 
sensitive to all these i nputs - as it occur s for the disk white dwarf 
luminosity function (llsern et al.ll2008h - and consequently that 
unless we have a more accurate determination of the luminos¬ 
ity function at large bolometric magnitudes (low luminosities) 
there is no hope to extract all this information from the observed 
data. However, large space-borne surveys, like Ga ia, will pro¬ 
vide us with a large sample of halo white dwarfs (iTorres et al.l 
l2005h . and hopefully all this wealth of information will be ex¬ 
tracted in a near future. Nonetheless, the lack of sensitivity of 
the hot branch of the luminosity function of halo white dwarfs to 
all these inputs can be interpreted in a positive way, since it al¬ 
lows us to obtain a robust statistical measure of the cooling rate 
of white dwarfs at low metallicities, and high luminosities. 
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